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Abstract 

We discuss temporal dark solitons in confined nonuniform Bose condensates. As a kind of local- 
ized high excitations, these solitons can be viewed as macroscopic quasiparticles, having relative 
motion to the background condensate. We get analytic expression for one dark soliton under slowly 
varying approximation and discuss its special propagation properties in nonuniform condensate, 
then we numerically prove that this approximation is reasonable and this kind of solitons exhibit 
their propagation properties in the nonuniform condensate. Finally, we simulate the generation of 
dark-soliton-like pulses in the condensate, and indicate that the excitation experiment, done by 
W. Ketterle's group can also be interpreted in terms of temporal dark soliton. 

PACS numbers: 03.75.Fi, 42.65.Tg, 42.81.Dp 
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Bose Einstein condensation(BEC) has been observed in dilute alkali atomic vapors JT], 
|2], ^]. Recently, further experiments have also demonstrated that the condensate can be 
perfectly described by the Gross-Pitaevskii equation ||, |I| §|, 0. For a ground-state conden- 
sate, the kinetic energy in this equation is usually negligible as compared with the s-wave 
scattering interaction energy that is the cubic nonlinear term. This indicates that the scat- 
tering term dominates main properties of the ground-state condensate, although it is a 
weakly interacting system. But for an excited condensate in which macroscopic number of 
atoms are collectively excited from ground state into high modes, the kinetic energy may 
become comparable to the scattering interaction energy. It could be noticed that the Gross- 
Pitaevskii equation, except for the confining potential, takes the form of the cubic nonlinear 
Schrodinger equation (CNLSE) which is well-known for the existence of soliton solutions ||. 
Naturally, one will ask whether solitons can exist in a condensate or not. In the previous 
works, A. Mysyrowicz et al. have already found the soliton-like propagation of the condensed 
excitons in CU2O crystal @, which obeys a similar Gross-Pitaevskii equation predicted by 
E. Hanamura | 10| . Weiping Zhang et al., G. Lenz et al., and P. A. Ruprecht et al. have 
proved the existence of bright soliton in coherent atomic waves [11, and condensate 
respectively. S. A. Morgan et al. have analyzed the solitary wave, synchronously moving 
with the background condensate The purpose of this paper is to discuss temporal dark 
solitons, as a kind of macroscopic quasiparticles, which have relative motion to the back- 
ground condensate, and it is different from the work in reference |T!| . This kind of solitons is 
actually localized high-mode excitations in the condensate consisting of atoms with positive 
scattering length, such as 87 Rb or 23 Na. First, we will consider a single temporal dark soliton 
in an idealized uniform condensate. Second, we will consider it in a nonuniform condensate 
and discuss its specific propagation features in this case. Third, we will numerically simu- 
late the propagation and collision of this kind of solitons in nonuniform condensate in order 
to prove that the dark solitons can exist in the nonuniform condensate steadily, exhibiting 
those propagation features, and the approximation adopted in the analytical solution is rea- 
sonable. Finally, we will also simulate the generation of dark-solition-like pulses by arbitrary 
perturbations, using the dark solitons' nonlinear superposition properties, and indicate that 
the excitation experiment, done by W. Ketterle's group H], can also be interpreted in terms 
of temporal dark soliton. 

For a highly excited condensate in which macroscopic number of atoms are excited into 
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high excited modes, the creation and destruction operators of these modes can be considered 



to be able to commute with each other |lo[ . Therefore, in an idealized uniform condensate, 
its atomic field operator can be treated as scalar, called macroscopic wave function. The 
macroscopic wave function can be described by the following Gross-Pitaevskii equation |R| . 



ih— ±±± = -— V 2 ^ (r, t) + Kxc* (r, t) + U | * (r, t)\ 2 * (r, t) (1) 

where \I/ (r, i) is the macroscopic wave function of the condensate, including not only the 
ground sate but also the macroscopically populated high excited modes. V exc is an external 
flat potential. Uq = 47th 2 a sc /m is a scattering constant, where a sc is the s-wave scattering 
length and m is the atomic mass. Eq. ([]]) is a CNLSE that possesses soliton solutions. 
When Uq > for positive scattering length, the solutions should be nonlinear superposition 
of dark solitons, according to the soliton theorems || . With the inverse scattering transform 
method M, we can solve this equation and get one-soliton solution, 



$ (r, t) = ^o 1 + ^*- iu *? e 21 ' e -<^c+t/ |*o| 2 >M ( 2 ) 



where \l/o is the background amplitude of ^ (r, t), fi s and v s are real constants, satisfying 
fi 2 + u 2 = 1, and /i s is called the eigenvalue of this dark soliton. In this equation, 



U Q \^ \ 2 m/h 2 k s ■ (r - r ) + fx s U \^ \ 2 t/h 



(3) 



where k s is a unit vector in the propagation direction of the soliton, ro is the center coordinate 
of the soliton. Therefore, the atomic number density distribution is 

|^(r,t)| 2 = \V \ 2 (l-v 2 s sech 2 T) (4) 

where v 2 is called the darkness of the dark soliton. Additionally, it is worthy to be em- 
phasized that dark soliton possesses no threshold, so it can be stimulated easily. According 



to the theorems of dark soliton [[HJ, an arbitrary perturbation of the wave function can 
be described as a nonlinear superposition of solitons, which indicates that soliton has its 
universality in the space described by the CNLSE ([I]). From expression (Q), we can find 
that the soliton is a localized function, because when z — > ±oo , its dark density decreases 
to zero quickly. Thus we can derive its full width at half peak darkness, 

arccoshV2 

A S (FWHM) = j = (5) 

v„\ 7ra sr l^nl 2 



This expression has a similar form with the healing length of the Bose condensate p8[ . 
except for the characteristic constant of the soliton. It is known that a two-dimensional 
vortex core in the Bose condensate has a similar size to the healing length. Consequently, 
their common features tell us that the dark soliton could be viewed as a kind of one- 
dimensional vortex core. Although there is no circular current around this kind of core 
due to its noncircularly-connected feature, the both sides of the soliton are connected by 
its bottom current density, which we will give more discussion later. Similar to a two- 
dimensional vortex, this one- dimensional vortex has a " vortex line" , or more exactly, it may 
be called vortex plane, which is the center plane of the soliton perpendicular to the vector k s . 
Corresponding to the circulation around a vortex line in two-dimensional case, there is also a 
similar quantity for the one-dimensional vortex core, which is just the multiply of the width 
and the velocity of the soliton and is directly determined by its characteristic constant. 
However, for a dark soliton, its characteristic constant is usually not a quantized, but a 
continuously valued number, this feature indicates that the one-dimensional vortex is not 
quantized for an uniform Bose condensate, and therefore the corresponding one-dimensional 
flow should exhibit classical fluidity instead of superfluidity. But, this does not mean that it 
won't become quantized for the externally confined nonuniform Bose condensate, because in 
considering the solitonic nonlinear dynamics at the vicinity of the confined boundaries, this 
problem always becomes very complicated and couldn't be solved analytically at present. 
With Eq. (^), we can calculate that when the atomic number density of the condensate 

|*o| is 1014 

cm 3 , the width of a black {y s = 1) soliton is 0.95/im for condensate of 23 Na. 
Obviously, it can be much smaller than the real size of a nonuniform condensate, for example, 
17/im in radial direction and 300/xm in axial direction, realized by W. Ketterle's group @. 
So the previous idealized uniform condensate is reasonable for very dark solitons. 

Subsequently, we consider a more real and special case in which the condensate, con- 
fined in external potential, is nonuniform, and it has a shape of cigar. Similar to the above 
assumption, a macroscopic number of atoms are still assumed to be excited, and the com- 
mutation relations of the atomic field operators are still valid. Thus, the Gross-Pitaevskii 
equation can be expressed as 

ih— ^ = -^V 2 * (r, t) + Kxt (r) * (r, t) + U |* (r, t) | 2 * (r, t) (6) 

where V cxt (r) = V r (R) + V a (z) is the external confining potential, which can be divided into 



two parts, the radial V r (R) and the axial V a (z) in cylindrical coordinates. Additionally, 
We assume the wave function \1/ (r, t) can be divided into the background ground-state 
condensate $ (r) and a soliton function j3 (z, t) propagating along z axis, 

* (r, t) = f3 {z, t) $ (r) e-^* (7) 

where \i c is the chemical potential of the stationary ground-state of the condensate. Then 
Eq. (||) is transformed into 



ih df3{z,t) h 2 d 2 p (z, t) h 2 dp {z, t) din® (r 



$(V) 



p(z,t) 



dt 2m dz 2 m dz dz 
+ [-U |/3 (z, t)®(r)\ 2 p (z, t) + U \$ (r) | 2 (3 (z, t)] ® (r 

u c ® (r) - ^-V 2 ® (r) + V ext (r) ® (r) + C/ |$ (r)| 2 ® (r) 
2m 



(8) 



Because ® (r) is the stationary ground-state wave function of the condensate, it must satisfy 



Vc® (r) 



2m 



V 2 ® (r) + V ext (r) ® (r) + U \® (r)| 2 ® (r) 



(9) 



Consequently, Eq. (RJ) is reduced to 



itj dP(z,t) + h 2 d 2 f3 (z, t) + h 2 dp (z, t) din® (r) 



dt 



2m dz 2 



m dz 



dz 



-U \P (z, t)®{r)\ 2 p (z, t) + U \® (r) | 2 P (z, t) = 



(10) 



Furthermore, we assume |$(r)| 2 is a slowly varying function of coordinate z, as compared 
with P (z, t) , satisfying 



d 2 p(z,t) 



dz 2 
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dp {z, t) din® (r) 



dz 



dz 



(11) 



Therefore, we can ignore the third term in Eq. (|10"D , and get 

dp(z,t) h 2 d 2 p(z,t) 



dt ~^ 2m dz 2 

i2 • 



U \P (z, t) ® (r)|" P (z, t) + U \® (r) \ z P (z, t) = 



;i2) 



Because |$(r)| is a slowly varying function in contrast to a soliton, when we solve this 
equation with inverse scattering transform method for a single soliton solution, we can take 
|$ (r)| 2 as a constant. Subsequently, we can routinely get one soliton expression, 

l + (^ 



P(z,t) 



iv>?e- 2T 



1 + e 



-2r 



(13) 



where 

T = u s ^U |$ (r)| 2 m/Vk s ■ (r - r ) + ti 8 U Q |$ (r)\ 2 t/h (14) 

fx s and v s still follow the previous definition, zq is the center coordinate of this soliton. The 
full width at half peak darkness of the soliton is 

arccoshy/2 

^s(FWHM) — > = [i-O) 
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which is different from Eq. (^) for |$ (r)| 2 being a slowly varying function of z. This means 
A 

s(fwhm) is to be increased when the soliton moves from the "top" to the "downhill" of the 
ground condensate wave function. From Eq. (|i~4j) we can derive the velocity of the soliton, 



v(D=- M |^k (16) 

m 

where k is the unit vector in the direction of z. v (r) varies with $ (r), which is consistent 
with the local speed of sound given by Bogoliubov |y|] and Lee, Huang, and Yang pO 



Additionally, Eq. ([16]) tells us that the velocity of this temporal dark soliton is also relevant 
to the eigenvalue, and it is, more explicitly, usually less than the absolute value of the 
corresponding sound speed in the nonuniform Bose condensate. However, this property can 
not be gotten in the density perturbation theory of hydrodynamics, such as ref. ||21|| . In 
considering the local density dependent property of the velocity, we can find that there 
is a slow feedback process between the velocity and the displacement of the soliton, and 



the temporal dark soliton behaves as an oscillator in the nonuniform Bose condensate [22 
Additionally, for those radial tightly confined Bose condensates, the velocity of the soliton 
varies rapidly with the density in the radial direction. This leads to very serious radial 
dispersion of the temporal dark soliton, to which a similar phenomenon has been observed 
in the experiments done by W. Ketterle's group |TJ]]. As it is known in nonlinear optics, the 



1 + 1 dimensional dark solitons are usually instable under their transverse perturbations, 
and always evolve into dark vortices f23fl . The analogy between this solitonic excitation in 
Bose condensate and the optical dark solitons enlighten us that the temporal dark soliton 
must be instable under the stretch of transverse dispersion, and it may also evolve into 
vortices. Due to the cylindrical symmetry of the ground state condensate, it is more likely 
to evolve into vortex rings [HJ, because they can still preserve the cylindrical symmetry. The 



circulation of those stable vortex rings must be quantized, which is a direct consequence of 



the coherent property of the Bose condensate, relevant to the nature of the superfluidity of 
the Bose condensate. Our further work on this topic is still in progress. Assuming that the 
radial variation could be neglected for some cases, we will find the product 

A s (fwhm) • v (r) = —1arccosh\ph^- (17) 

should be a constant for a darkness-fixed soliton. As we have described above, analogous to 
the circulation around the vortex line of a two-dimensional vortex in the fluid, this quantity 
may also be called a one- dimensional circulation, which characterizes the properties of the 
one-dimensional dark soliton vortex core directly. But the more important thing is that 
it is very realistic for an experimental measurement of this quantity. Therefore, Eq. (|17D 
may offer us a useful method to judge whether a dark pulse in the condensate can be very 
precisely interpreted as a temporal dark soliton. From Eq. (^) we can also get the atomic 
current density carried by the soliton 

j = rq,* ( r t ) W (r, t) - V (r, t) W* (r, t)) = kfi s v 2 Un /2 m- 1/2 $ 3 (r) sech 2 T (18) 

2m 

which is pulsed. The width of this pulse is proportional to that of the soliton. It is worthy 
to be noticed that the direction of the current density pulse is opposite to the propagation 
direction of the dark soliton, however, the propagation direction and the velocity of the 
current density pulse are both same to those of the soliton. When v s — or v s — 1, the 
amplitude of this current density pulse reaches its minimum j min = 0. When v s = 2^3/9, it 
gets to its maximum 

W = ^U^m-^ (r) (19) 

The existence of this maximum shows that each distribution of the current density usu- 
ally corresponds to two possible soliton eigenvalues and therefore two one-dimensional dark 
soliton vortices. This property is different from the usual two- dimensional vortex case. 
Additionally, from Eq. (]18|) , we can find that the relationship of the current density and the 
soliton eigenvalue fi s is in the shape of S. In considering some nonlinear dynamical process, 
such as the solitonic excitation or the soliton collision with a boundary. This relationship 
may provide us an opportunity of finding bistable phenomena of the temporal dark soliton. 

In the above derivation, we have adopted slowly varying approximation (|TT|). To make 
sure this approximation is reasonable and temporal dark solitons can really exhibit their 
propagation properties in a confined nonuniform condensate, we numerically simulate the 
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propagation and collision of this kind of solitons. We assume the ground-state condensate 
is composed of N = 5 x 10 6 of 23 Na atoms with a sc = 27. 5A |[24|| , and the radial trapping 
frequency is u r = 2tt x 1800Hz, the axial trapping frequency is u a = 2ir x 18Hz. During the 
simulation, we adopt one-dimensional Thomas-Fermi approximation in the radial direction 
of the cylindrically symmetric condensate. Consequently, we can approximate Eq. ([]) as 
one-dimensional Gross-Pit aevskii equation by integration in the radial cross section of the 
condensate, 

*ft g- t + ^ctpa{z,t) = - — — + V a (z)ij a (z,t) + -Uo\lp a (z } t)\ 1p a {z,t) (20) 

where ip a (z, t) is the axial wave function of the condensate, distributed along the axial direc- 
tion z, and V a (z) = muj 2 z 2 /2. Although the one-dimensional Thomas- Fermi approximation 
is still not strict, it can still keep the axial wave function ip a (z,t) without slowly varying 
approximation ([0]) in the direction of z, and therefore Eq. (|20D is sufficient to give the ap- 
proximation a test. We numerically solve Eq. fl20|) by split-Fourier transform method [p5 |. As 
an example, two solitons' collision process is illustrated in Fig.l. The two soliton, v s \ = 0.5, 
v S 2 = 0.5, set initially on the "downhill" of the condensate wave function with distance 
6.2/xm, shown in Fig.l. (a). It should be noted that both the wave function and the atomic 
number density, shown in all figures in this paper, have already been transformed into di- 
mensionless forms with respect to their ground-state center values at z = 0.0/im. These 
two solitons propagate and collide with each other, and then they separate and propagate 
independently, as shown in Fig.l.(b). Before and after the collision, they can conserve their 
shape steadily. From the contour graph of this collision, as shown in Fig.l.(c), the widths 
of the solitons vary with their position, which indicates the locality of their widths. The 
velocities of the two solitons are all about v = 15fim/ms. There is no damage occurred to 
both of these solitons, so the collision must be elastic. More numerical experiments show 
that so long as the numerical precision is enough, the solitons can propagate steadily as far 
as one can compute. Therefore, we can draw a conclusion that the slowly varying approxi- 
mation (|TT| ) is reasonable, and the solitons indeed have the propagation properties, indicated 
above in nonuniform condensate confined by external potential. Additionally, as we have 
seen, the properties of these temporal dark solitons in the condensate are very similar to 
those of conventional particles, so we can consider them as a kind of quasiparticles. But they 
are different from phonons, because every one of them is macroscopic whereas a phonon is 
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usually microscopic, and furthermore every soliton must be a coherent structure as described 
in Eq. QT|). 

According to the dark soliton theorems ||, dark solitons have nonlinear superposition 
properties. We may use these properties to generate temporal dark solitons in nonumiform 
condensate. To illustrate this possibility, we numerically simulate the generation of soliton- 
like pulses by arbitrary perturbations in the nonuniform condensate. As an example, here 
we give two reversed Gaussian function perturbations in the condensate, shown in Fig.2.(a). 
The darkness of these two perturbations are 0.44 and 0.64 respectively, and the full widths 
of them are 2.47/xm and 0.99/im respectively, and their distance is 10/im. As shown in 
Fig.2.(b), each of these two dark perturbations splits quickly into more pairs of dark pulses, 
propagating with different velocity depending on its own darkness. The collision of the two 
out-splitted dark pulses must be elastic because the two pulses can still keep even symme- 
try with their own twins. Obviously, these properties are very similar to those of solitons. 
From Fig.2.(c), which is the contour graph of Fig. 2. (a), we can find that all the velocities 
of these out-splitted soliton-like pulses are around V = 20fim/ms. Although we can't make 
sure that these out-splitted pulses won't split anymore ( in fact, the pair, splited by the 
left one, have already begun to split again into more pairs of pulses), with the above elastic 
collision evidence, we can at least draw a conclusion that the initial dark pulses can be 
approximated as multi-solitons, each of which is composed of a number of solitons. The 
recent experiment, done by W. Ketterle's group has shown a similar phenomenon as 
we have simulated here. In their experiment, one dark perturbation splits into two, and 
then these two out-splitted pulses continue to propagate with spreading widths. Although 
there is no description about observation of further splitting of the out-split dark pulse, the 
out-split dark pulse must be composed of more than one solitons. The spreading of the 
pulses can be interpreted as the radial dispersion as we have described above, but one can't 
exclude some other possible reasons, such as further splitting of the multi-solitons and the 
widening with the decrement of the background condensate as described by Eq. (|T5|). And 
it can't be usually simply explained only as the last one, except that one can, by chance, 
get twin strict dark solitons, every one of which must have the form as we have described 
analytically in Eq. flT3|). However, further splitting may not be observable, because the size 
of the condensate is limited, the velocities differences are all to decrease to zero due to the 
Eq. (|I~6|), and they may therefore have no enough distance to propagate for further splitting. 
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To test the soliton-like properties of the dark pulses in the present experiment, according to 
our numerical simulation, we suggest that one can generate two dark pulses propagating in 
opposite direction, and observe their collision. One can also test these properties by measur- 
ing the relationship between the velocity of every steadily propagating out-split dark pulse 



tion of the experiment would need more study in details and more realistic simulation which 
can't be included in this short paper. But, as we have briefly discussed above, the temporal 
dark soliton theorems are promising tools to interpret the experiment done by W. Ketterle's 



In conclusion, we have derived the temporal dark soliton in nonuniform condensate under 
slowly varying approximation. Viewed as the one-dimensional vortex cores, these temporal 
solitons may provide deeper insight into the superfluidity in the cigar-shaped Bose-Einstein 
condensate. Then we have numerically simulated the propagation and collision of this kind 
of solitons, and proved that the slowly varying approximation is reasonable and the solitons 
can exhibit the special propagation properties in the nonuniform condensate. We have 
finally simulated the generation of soliton-like pulses by two arbitrary dark pulses in the 
nonuniform condensate, and with this simulation, we have given the experimental results of 
W. Ketterle's group a qualitative interpretation in terms of temporal dark soliton. 
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Fig. 1 a "The Propagation of Sound Dark Solitons in Bose-Einstein Condensates" 
T. Hong, et al. 




Fig. 1 b "The Propagation of Sound Dark Solitons in Bose-Einstein Condensates" 
T. Hong, et al. 
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Fig. 1 c "The Propagation of Sound Dark Solitons in Bose-Einstein Condensates" 
T. Hong, et al. 



FIG. 1: The collision of two temporal dark solitons in the condensate, (a) The dimensionless 
axial wave function of the ground- state condensate is initially excited with two dark solitons. (b) 
The two temporal dark solitons propagate steadily and collide with each other elastically. (c) The 
contour graph of the two colliding solitons. Both the widths of them vary with their position on 
the condensate. 
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Fig. 2 a "The Propagation of Sound Dark Solitons in Bose-Einstein Condensates" 
T. Hong, et al. 
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Fig. 2 b "The Propagation of Sound Dark Solitons in Bose-Einstein Condensates" 

T. Hong, et al. 
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Fig. 2 c "The Propagation of Sound Dark Solitons in Bose-Einstein Condensates" 
T. Hong, et al. 

FIG. 2: The excitation of dark temporal dark soliton in a ground- state condensate. 

(a) The dimensionless axial wave function of the ground- state condensate is initially excited with 
two reversed Gaussian shape dark pulses. 

(b) Each of two dark pulses splits into pairs of dark pulses, and the split pulses collide with each 
other elastically. This indicates that the initial dark pulses can be approximated as nonlinear 
superposition of temporal dark solitons. 

(c) The contour graph of the two dark pulses. Some of the split dark pulses propagate steadily 
and some of them are splitting into more dark pulses. 
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